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Abstract 

In Moyal's formulation of quantum mechanics, a quantum spin s is described 
in terms of continuous symbols, i.e. by smooth functions on a two-dimensional 
sphere. Such prescriptions to associate operators with Wigner functions, P- or 
Q-symbols, are conveniently expressed in terms of operator kernels satisfying the 
Stratonovich-Weyl postulates. In analogy to this approach, a discrete Moyal for- 
malism is defined on the basis of a modified set of postulates. It is shown that 
appropriately modified postulates single out a well-defined set of kernels which 
give rise to discrete symbols. Now operators are represented by functions taking 
values on (2s -|- 1)^ points of the sphere. The discrete symbols contain no redun- 
dant information, contrary to the continuous ones. The properties of the resulting 
discrete Moyal formalism for a quantum spin are worked out in detail and com- 
pared to the continuous formalism, and it is illustrated by the example of a spin 
1/2. 

1 Introduction 

The idea to represent quantum mechanics of a particle in phase space F goes back to 
Wigner He established a one-to-one correspondence between a quantum state lip) 
in the particle Hilbert-space Ti. and a real function 

2 f 

W^{p,q) = —j dxilj*{(l + x)ilj{q — x) exp[2ipx/h] . (1) 

Its properties suggest an interpretation as a quasi-prohahilitj in phase space, the only 
'drawback' being due to the negative values it may take. A more general framework for 
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phase-space representations of quantum states as well as operators A is given by the 
relation 

W^iq,p)= Tr 

with an operator kernel 0, ^ 

A{q, p) = 2D{q, p) fi D\q, p) , (g, p) G T . (3) 

Here H : {q,p) i—Q,—p) is the unitary, involutive parity operator while D{q,p) 
describes translations in phase space 0. If the operator A is chosen as the density 
matrix of a pure state, A = = \ip){ip\, Eq. (|D reduces to (0). The kernel A{q,p) 
can be derived from a set of conditions which a phase-space representation is required 
to satisfy (cf. below). It is intimately related to the behaviour of a function WA{q,p) 
under translations mapping the phase space F onto itself. The map from operators 
to functions {A Wa) has an important feature: its inverse, mapping functions to 
operators (Wa A), is mediated by the same kernel — in other words, the kernel 
A{q,p) is self-dual. 

For a quantum spin, the symbol associated with an operator is a continuous function 
defined on a sphere iS^, which is the phase space of classical spin. Now, instead of trans- 
lations in planar phase space, it is the group SU (2) of rotations which plays a dominant 
role when the Moyal formalism is set up. As for a particle, the set of Stratonovich-Weyl 
postulates characterizes the symbols in an elegant way. For clarity, the postulates 
are now displayed in their familiar form for the continuous symbols: 

(5*0) linearity : A Wa is a linear one to one map, 

(SI) reality: W^^t (n) = W^l(n) , 



A{q,p)A 



(2) 



(S2) standardisation: 2s + 1 /" , , , m !7^ 

{S3) traciality : WA{n)WB{n)dn = Tr [AB] , 



(54) covariance : Wg.A = W%, gE SU{2) . 

It is natural to have a linear relation between operators and symbols (SO), while (SI) 
implies that hermitean operators are represented by real functions. The third condition 

(52) maps the identity operator to the constant function on phase space, and traciality 

(53) ensures that the correspondence between operators and symbols is invertible. The 
covariant transformation of the symbols with respect to rotations g G SU{2) effectively 
introduces phase-space points as arguments of the symbols. The continuous Moyal 
representation for a spin 0, || compatible with these conditions can be based on a 
self-dual kernel A(n) (cf. Sect. 2) in analogy to (Q). 
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In order to have a consistent and full-fledged classical formalism it is necessary to 
introduce a product between symbols which keeps track of the non-commutativity of 
the underlying operators. This Moyal product p[, or twisted product, for two operators 
A and B expresses the Wab of the operator product AB in terms the symbols Wa and 

WABin) = WA(n)*WB(n)= [ [ L(n, m, k)W^^(m)M^B(k)dmdk , (4) 

with a function L(n, m, k) of three arguments given explicitly in for example. The 
* product is known to be associative. 

Other continuous representations for a spin do exist, such as the Berezin symbols 
of spin operators which are the analog of the P- and Q-symbols for a particle. 
Instead of a single self-dual kernel, the Berezin symbols require however, a pair of two 
different kernels, dual to each other: one of the kernels maps operators to functions 
while its dual is needed for the inverse procedure. It will become clear later on that the 
self-dual and the dual approach correspond to deflning orthogonal and non-orthogonal 
bases, respectively, in the vector space As of operators acting on the Hilbert space of a 
spin s. When slightly modifying the postulates of Stratonovich, they are also compatible 
with kernels which are not self-dual. 

A common feature of these representations is the redundancy of the continuous sym- 
bols. When represented by a (2s + 1) x (2s + 1) matrix, a hermitean operator is flxed 
by the values of (2s + 1)^ real parameters. Consequently, the values of the symbols, 
continuous functions on the sphere, cannot all be independent — in other words, the in- 
formation contained in a symbol is redundant. The discrete version of P- and Q-symbols 
for a spin s, introduced in |Tl| as a means to reconstruct the quantum state of a 
spin, allows one to characterize a spin operator A by using only the minimal number 
of parameters. In fact, a discrete symbol can be considered as living on a 'discretized 
sphere,' that is, as a function taking (real) values on a flnite set of points on the sphere 
only. Such a formalism will be called a discrete Moyal-type formalism. 

The purpose of the present paper is to develop the discrete Moyal formalism in 
analogy to the continuous one. In particular, the kernel and its dual deflning the dis- 
crete symbols will be derived from a set of appropriate Stratonovich-type postulates. 
Subsequently, the properties of these symbols are studied in detail. 

2 Continuous representations 

Continuous self-dual kernel: Wigner symbols 

The Stratonovich- Weyl correspondence for a spin s is a rule associating with each oper- 
ator A E As on a Hilbert space Hs a function Wa on the sphere S^, called its (Wigner-) 
symbol. Let us deflne it in analogy to Eq. @, by means of a universal operator kernel 
A(n), which can also be thought of as a fleld of operators on the sphere. Then, the flrst 
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requirement (SO) is already satisfied, and the postulates (SI) to (S4) turn into conditions 
on the kernel: 

(CI) reality : At(n) = A(n) , 



(C2) standardisation : 2s + I 



dnA(n) = I . 



4tt Js^ 

iC3) traciality: rfn Tr [A(n)A(m)] A(n) = A(m) , 

(C4) covariance : A(^ ' n) = f/^ A(n) f/t , ge SU{2) . 

where the matrices Ug are a unitary (2s + 1) -dimensional irreducible representations of 
the group SU (2). 

The existence of a kernel A(n) satisfying (Cl-4) has been proven in by explicit 
construction. The derivation in starts by expanding the kernel in a basis associated 
with the eigenstates of the operator § ■ n^, 

s 

A(n) = Zw(n)|m, n^)(m',n^| , (5) 

•m,m' =—s 

with unknown coefficients ^^^'(n). It follows from (Cl-4) that one must have 

(n) = ^^ye;V27TT/ / 
mm\ ) 2s + lf^^ \ m m' -m 



^/ y>^,m'-m(n), (6) 



where Eq = 1 and ei = ±1 , / = 1, ... ,2s, and the definition of Clebsch-Gordan coeffi- 
cient given in is used. Consequently, there are 2^* different kernels which define a 
Stratonovich-Weyl correspondence rule. 

A new and simple derivation of the kernel A(n), independent of the argument given 
in 0, is presented now which has two important advantages. On the one hand, it will 
provide a form of the kernel similar to that one of a particle (^, which is interesting 
from a conceptual point of view. On the other hand, it will be possible to transfer this 
approach to a large extent to the case of the discrete Moyal formalism. 

Expand the kernel in the eigenbasis of the operator s ■ n, 

s 

A(n) = Amm/(n)|m, n)(m', n| , (7) 

m,m'=—s 

where the expansion coefficients A^m' are unknown so far. According to the reality 
condition (CI) they must satisfy Amm'i'^) = Aj^,„(n). In a first step, the numbers 
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Amm'(ii) are shown not to depend on the label n. Consider the transformation of A(n) 
under a rotation g. According to (C4) one must have 

s 

A,w(fl' ■ n)|m,5( ■ n)(m',5( ■ n| = 

s 

UgA{ji)Ul= J2 ^mm'in)\m,Rgn){m',Rgn\, (8) 



m,m——s 



m,m' =—s 



where Ug\m, n) = |m, -R^n) = |m, g ■ n) with a rotation matrix Rg G 5*0(3) representing 
g G SU{2) in IR^. Consequently, one must have 

Aw(5( • n) = Aw(n) , (9) 

which is only possible if A^m' does not depend on n. Consider next a rotation (^(n) about 
the axis n by an angle G [0, 27r), represented by the unitary f/g(n) = exp(m ■ sip). The 
left-hand-side of (C4) is invariant under this transformation while the right-hand-side 
transforms: 

s 

A(i?g(n)n) = A(n) = ^ Aw exp[z(m - m')v9]|m, n)(m', n| , (10) 

m,m'=— s 

which is possible only if 

Amm' = A{m)6mm' ■ (H) 

Therefore, covariance of the kernel under elements of SU (2) requires it to be diagonal 
in the basis associated with the direction n, 

s 

A(n) = ^ A(m) |m, n) (m, n| . (12) 

m=—s 

Next, the condition of traciality will be exploited. Upon rewriting (C3) in the form 

/ rfn(5,(m,n)A(n) = A(m), (13) 

the function (5s(m, n) = (2s + 1) Tr [A(m)A(n)]/(47r) is seen to be the reproducing 
kernel for a certain subset of (2s -|- 1)^ functions on the sphere [|12|, 0. In other words. 



5s(m, n) acts in this space as a delta-function with respect to integration over 5^, and 
for spin s, it reads explicitly 



2s I 2s 



21 4- 1 

<5.(m, n) = E 5: YUni)YC^{n) = ^ -^Pi{m ■ n) . (14) 

1=0 m=-l 1=0 
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Here the addition theorem for spherical harmonics Yim{n), I = 0, . . . ,2s, —I < m < I, 
has been used to express the sum over m in terms of Legendre polynomials Pi{x), —1 < 
X <1. Upon choosing m = and with ■ n = cos 6', the condition (C3) becomes 



Tr 



AfnJAfnl 



, , P;(cos( 



1=0 



2s + 1 



(15) 



Use now the expansion (^) of the kernel on the left-hand-side as well as the identity 



2s 



|(m,n2|m', n)|^ = ^ 



z=o 



/ s I 




1 s I 


s 


\ m 


m 


\ m' 


m' 



Piicose), (16) 



leading to 



2s 



Tr 



A(n.)A(n) = E E ^1 



m] 



1=0 \m=—s 



S I 

m 



s 

m 



21 + 1 
2s + 1 



Pi{cos9). (17) 



Compare now the coefficients of the Legendre polynomials Pz(cos^^) with those in Eq. 
([151) . This leads to (2s + 1) conditions 



E A( 



m] 



s I 
m 



s 

m 



ei = ±l, 1 = 0,..., 2s. 



These equations can be solved for A(m) by means of an orthogonality relation for 



Clebsch-Gordan coefficients 13 



A m = > El ( 

^ ' ^ '2s + 1 \ m 



1=0 



s 

m 



(19) 



Thus, the self-dual kernel for the continuous Moyal formalism is given by 



A(n)= E E^/ 



2/ + 1 / s / 



m=—s 1=0 



2s + 1 \ m 



m 



\m, n) {m, n| 



(20) 



Out of these 2^*"'"^ distinct solutions only 2"^^ are compatible with the condition of stan- 
dardization (C2) which has not been used until now. This condition imposes 



E A( 



m] 



(21) 



being satisfied if and only if Eq = +1. 

The set of solutions (pOD coincides indeed with those found in 0. The easiest way 
to see this is to calculate the matrix elements of the kernel A(n) in ( pOD with respect 
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to the standard basis |m, n^). One reproduces the coefficients of the expansion (j^): 
(m, n^|A(n)|m',n2) = Zw(n). 

The expansion ( pOD is interesting from a conceptual point of view. It allows one to 
interpret physically the kernel A(n) in analogy with the kernel for a particle given in 
(0) by writing 

A(n) = [/„A(n,)f/t, (22) 
where represents a rotation which maps the vector on n. Imagine now to contract 



T2| the group SU (2) to the Heisenberg-Weyl group. It is known this procedure turns 
rotations Ug into translations D{q,p) . As shown in |T^, the operator ^{11^) contracts 



in the following way, 

A(n,)^2n, (23) 

if 61 = +1,1 = 0,. . . ,2s. Therefore, the operator A(n2) plays the role of parity for a 
spin which is by no means immediately obvious when looking at it. 

Finally, we would like to point out that the integral kernel L, defining the * product 
of symbols (§), has a simple expression in terms of Wigner kernels: 

L(n,m,k) = C^^^) Tr [A(n)A(m)A(k)l . (24) 
V 47r / L J 



Continuous dual kernels: Berezin symbols 

Wigner symbols of spin operators are calculated by means of a kernel which is its own 
dual. Other phase-space representations are known which do not exhibit this 'symmetry' 
between an operator and its symbol. P- and Q-symbols for a particle are familiar 
examples which have their analog in the 'Berezin' symbols for a spin. It will be shown 
now that these symbolic representations also have a simple description in terms of kernels 
satisfying a modified set of Stratonovich-Weyl postulates. The conditions (Cl-4) must 
be relaxed slightly in order to allow for a pair of dual kernels. 

The required generalization is easily understood in terms of linear algebra. The 
ensemble of all operators, that is, the self-dual kernel is nothing but a an (overcomplete) 
set of vectors spanning the linear space As of operators on the Hilbert space of the spin 
s. As the traciality (C3) indicates, this family of vectors is 'orthogonal' with respect to 
integration over the sphere as a scalar product. Each operator A can be written as a 
linear combination of the elements of the kernel with its Wigner symbol as expansion 
coefficients. More precisely, the expansion coefficients 1^^(11) with respect to the basis 
A(n) are given by the 'scalar product' of A with the same basis vector as shown, for 
example, in Eq. (0). The essential point now is, that there are also non-orthogonalhases 
of the same space. Given a non-orthogonal basis, denoted by An, its dual basis A" 
is uniquely determined through the scalar product. Furthermore, the dual basis also 
spans the original space which implies that now there will be two different expansions 
of one operator A defining a symbol A^ and its dual A". Consequently, both kernels 
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and symbols now come in pairs. The familiar P- and Q-symbols — or Berezin symbols 
— will turn out to be related in this precisely way. 

Non-orthogonal bases are allowed in the present framework if, first of all, traciality 
(C3) is relaxed to 

2s + 1 



(C"3) traciality : 



dn Tr 



An = A, 



(25) 



The kernel and its dual are both real in analogy to (CI). Explicitly, the symbols and 
their duals are given by 



Tr 



AAr 



A"" 



Tr 



AA' 



(26) 



Furthermore, one is free to normalize one of the kernels. An, say, in analogy to (C2), 
and one requires it to transform covariantly (C4). 

It is possible as before to derive the explicit form of the kernels by a reasoning in 
analogy to above. The general ansatz for both An and A" in the basis referring to the 
axis n as in (|^) is again reduced to diagonal form by exploiting their behaviour under 
rotations: 



Am|^, n)(m, n| , A" = ^ A^lm, n)(m, n| 



(27) 



with two sets of numbers A^ and A™, which do not depend on n. It is necessary that 
the trace of these two operators with labels m = and n, say, equals the reproducing 
kernel with respect to integration over the sphere, that is, instead of (|l^) one needs to 
have ^ 

Tr [An.A"l = ^^P/(cos^) , (28) 



=0 



2s + 1 



where cos^ = ■ n. This leads to the conditions 



s 

m 



s 
m 



s I 
m 



s 

m 



(29) 



with / = 0, . 
numbers 7/: 



2s. The ensemble of solutions is parameterized by (2s + 1) non-zero real 



E A 



/ s / 
""^ m 



s 

m 



E 



s / 
m 



s 

m 



7« 



Solving for the expansion coefficients, one obtains 



Ajn 

A™ 



^ 21 + 1 



I 



2s + 1 \ m 



^ _i2/ + l 
1=0 



s 

m 



2s 



s I 
m 



s 

m 



(30) 

(31) 
(32) 
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As in the self-dual case, the standardisation implies that 70 = +1. This class of solutions 
for kernels which are not their own dual has been obtained in by an entirely different 
approach. Self-dual kernels are a small subset: they require = A*", which is 7/ = 7^"^ 
or 7; = ±1 in agreement with (pOf ). Each set of numbers 7; defines a consistent phase- 
space representation of a quantum mechanical spin. 



Consider the particular case A^ = 5ms resulting from 



71 



s I 
s 



The associated kernels read 



A, 



A" 



|s, n)(s, n| = |n)(n| 



2s 



E E 

m=—s 1=0 



1 s I 


s ' 


)\ 


' s I 


s 


\ s 






m 


m 



|m, n)(m, n| ., 



(33) 

(34) 
(35) 



This choice has the advantage that one of the two symbols, reducing to the expectation 
value of an operator A in coherent states 0, is particularly simple. It turns out to be 
just the Q-symbol, Qyi(n) = (n|A|n), that is, its expectation in a spin-coherent state. 
At the same time, one falls back on a familiar expression for the dual symbol which turns 
out to be the P-symbol for A, defined by an expansion in terms of a linear combination 
of operators projecting on coherent states. 



A 



(2.3 + 1] 



PA(n)|n)(n|(in. 



In the present notation one simply has in view of (120) that 

Q^(n) =A^ = Tr [iA„] , P^(n) = A^ = Tr [AA^ 
so that (|36f) reads 



A 



(23 + 1) 



477 

It is obvious now that one has (cf. 0) 



rfn Tr [IA"] a. 



Tr [AB] 



(23 + 1) 

477 



dnPA{n)QB{ri) 



(2g+ 1) 

477 



(36) 

(37) 
(38) 

(39) 



Finally, it is interesting to calculate the Q- and P-symbols of the self-dual kernel A(n) 
as well as the pair of dual kernels A" and An using the short-hand 



Yzm(n, m) 



(40) 
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so that the reproducing kernel is given by 

5,(n,m) =^Y;^(n,m). (41) 

ml 





TV [ • A„.] 


Tr [ • A'"] 


An 




5,(n,m) 


A" 


5s{n, m) 




1 S I 

y m 


4^ Yim(n,m) 


A(n) 


^-'(m ^)Wn,m) 


Y^ral { 


/ s I 
y m 


4^ Y/m(n,m) 



Note that the entries of last row, the Q- and P-symbols of the self-dual kernel A(n), do 
simultaneously provide the Wigner symbols of the dual kernels Am and A™. 

3 Discrete Moyal-type representations 

A particular feature of the kernels discussed so far is their redundancy: the linear 
space of hermitean operators for a spin s has dimension (2s + 1)^ while the kernels 
consist of a continuously labeled set of basis vectors. In other words, there are at most 
Ng = (2s + 1)^ linearly independent operators among all A(n), n e S'^. In this section 
discrete kernels will be introduced, denoted by Aj,, v — 1, . . . ,Ns. No linear relations 
must exist between the operators A^ which constitute the kernel, that is, they are a 
basis of As in the strict sense. It is natural to expect that a subset of precisely A^^ 
operators A(n,^), u = 1, . . . , Ng will give rise to a discrete kernel. Therefore, evaluating 
a continuous symbol of an operator A at Ng points rij, of the sphere S^, provides a 
promising candidate for a discrete symbol, i.e. the set A^, = A-n^, v — 1, . . . , A^^. For 
brevity, A^^ points on are called a constellation. 

As before, one might expect orthogonal and non-orthogonal kernels to exist. It turns 
out, however, that an appropriately modified set of Stratonovich-Weyl postulates cover- 
ing discrete kernels does not allow for orthogonal ones. Therefore we start immediately 
by deriving the discrete non-orthogonal kernels coming as before in combination with a 
dual. 
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Discrete dual kernels 



By analogy with the continuous representation of the preceding section, one modifies 
the Stratonovich-Weyl postulates in the following way (throughout the index v takes all 
the values from 1 to A^J: 



(DO) linearity : 
{D\) reality : 

(D2) standardisation 



(D3) traciality 



Ay is a linear map . 
At=A,, z/ = l,...,Ar,, 



1 



2s + 



1 .4:1 



2s + 



-y Tr 

1 



A, 



(D4) 



covariance 



A,., = t/,A,t/t g^SU(2) 



Let us briefiy comment on these conditions. Linearity is automatically satisfied if discrete 
symbols are defined via kernels, that is. Ay = Tr [AA;^]. The second condition, reality, 
is obvious, and in (D2) the kernel A'" dual to A^, is standardized. The duality between 
Ay and A'^ is made precise by the condition of traciality since (D3) only holds if one 
has 

1 



Tr 



A.A^ 



Z/,/i 



(42) 



2s + 1 

which, upon considering the trace as a scalar product, is precisely the condition defining 
the dual of a given basis. As a matter of fact, if {A;^}, z/ = 1, . . . , A^^, is a basis, its 
unique dual is guaranteed to exist. Finally, covariance under rotations g G SU{2) must 
be reinterpreted carefully. Under a transformation g, a constellation associated with A^^ 
points on the sphere will, in general, be mapped to a different constellation. In other 
words, the image Ag.y = A{g-ny) is typically not one of the operators A^. Nevertheless, 
condition (D4) is not empty: for appropriately chosen rotations g^^ one can indeed map 
an operator defined at to another one associated with the point n^, say. In this 
case, the consequences for the coefficients of the operators Ay and A^ are identical to 
those obtained in the continuous case. Similarly, invariance of the operator A^, under a 
rotations about the axis iiy has the same impact as before. Thus the general ansatz for 
the discrete kernel (obtained from (0) by setting n n.y) is reduced by exploiting the 
postulates (Dl-4) to the form 



A, = A{ny) 



E 



Am\m,ny){m,ny 



V 



AT, 



(43) 



Therefore, the discrete kernel A^ can be thought of as a subset of A''^ operators A(ny), 
each one associated with a point \\.y of the sphere. 
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Let us mention an important difference between discrete and continuous kernels, Ky 
and A(n), wliicli arises in spite of tlieir formal similarity. Once the coefficients Am are 
fixed a continuous kernel is determined completely. Discrete kernels, however, come in 
a much wider variety since they depend, in addition, on the selected constellation of 
points on the sphere. The discrete kernel does not enjoy the SU{2) symmetry in the 
same way as does the continuous one. The discrete subgroups of SU{2) being limited 
in type, the continuous symmetry will usually not be turn into a discrete one. Note, 
further, that the elements of the dual kernel depend, in general, on all the points of the 
constellation: = A'^(ni..., n^r^). This is easily seen from (^2p since the variation of a 
single Ay will have an effect on all {A'^} in order to maintain orthogonality. 

The additional freedom of selecting specific constellations is connected to a subtle 
point: actually, not all constellations of A^^ points give rise to a basis in the space As- 
This remark is easily understood by considering as an example of a linear space. The 
(continuous) collection of all unit vectors in three-space clearly spans it while not every 
subset of three vectors is a basis — they might lie in a plane. By analogy, one must ensure 
that the operators Aj,, u = 1, . . . ,Ns, associated with a specific constellation, do indeed 
form a basis of As- The operators are indeed linearly independent if the determinant of 
their (positive definite and symmetric) Gram matrix G [O] satisfies 



det G > , 



Tr 



(44) 



a condition, which will be studied later in more detail. 

Suppose now that the A^^ operators in ( ^31) do form a basis. Then, the kernel 
dual to it, that is the set of operators A'', is determined by the condition (^2|) instead 
of Eq. (p!7D. Therefore, one cannot proceed as before to derive the conditions (p9|). In 
particular, it is no longer true that the elements of the dual kernel have an expansion 
analogous to (j43|). This follows immediately from the impossibility to satisfy (^) by 
an ansatz for A'^ of the form (|43|) : Eq. (^31) represent Ng conditions but a dual of the 
form depends only on (2s + 1) free parameters A™. Nevertheless, a dual kernel A'^ does 
exist and it is determined unambiguously — it simply cannot have the form (|43|) . [A 
also ori, there is no self-dual kernel associated with the Stratonovich-Weyl postulates 
(DO-4)). Consequently, one expands any (self-adjoint) operator A either in terms of a 
given kernel. 



A 



2s + 



A" 



Tr 



u=l 



or, equivalently, in terms of the dual kernel. 



A 



2s + 



Ay= Tr 



AA' 



AA„ 



(45) 



(46) 



The collection A = {Ai, . . . , A^J of real coefficients in (|4^) now is defined as the discrete 
phase-space symbol of the operator A, and A'^""' = {A^, . . . , A'^^) is the dual symbol. 
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The relation between the discrete symbol and its dual as well as between the pair 
of kernels is linear. It is easily implemented by means of the Gram matrix G and its 
inverse G~^, 

1 



(G-^) =0""' 



(2s +1)2 

The matrix G thus plays the role of a metric, 



Tr 



A^ = (2s + 1)$:G'^''A^, 
and the dual symbol is determined according to 

= {2s + ■ 

u'=l 



(47) 



(48) 



(49) 



The trace of two operators A and B is easily found to be expressible as a combination 
of a discrete symbol and a dual one, 



Tr 



AB 



Y,A''B, = Y,A,Br 



(50) 



v=l 



which is the discretized version of Eq. (|39|) . 

In order to have a discrete Moyal product, we seek to reproduce the multiplication 
of operators on the level of symbols. Using the definition of the symbols, it is straight- 
forward to see that 



{AB)x = Ax*Bx 
with the trilinear kernel 
in close analogy to Eq. (W). 



1 



Ns 



Lf = Tr 



A^A^A) 



(51) 



(52) 



Discrete P- and Q-symbols 

A particularly interesting set of symbols emerges if, for a given allowed constellation, 
only one of the coefficients in the expansion (|43|) is different from zero. Am = 6ms, say. 
Then, the kernel consists of A^^ operators projecting on coherent states, 

= |n,,)(ni,| . (53) 
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This is obviously the non- redundant counterpart of Eq. ( p^D implying that a self-adjoint 
operator A is determined by a symbol which consists of A^^ pure-state expectation values, 
the discrete Q-symbol, 

A, = Tr [AQ,] = (n,|i|n,) . (54) 

Let us point out that the introduction of discrete symbols has actually been triggered 
by the search for a simple method to reconstruct the density matrix of a spin through 
expectation values 0]. In fact, this problem is solved by Eq. (|54D in the most economic 
way. If A is chosen to be the density matrix p of a spin s, then the i^-th component of 
the Q-symbol equals the probability of measuring the eigenvalue s in the direction rij,, 

Psin^,) = {n^\p\n^) . (55) 

Knowledge of the Ng measurable probabilities Psi'n.u) thus amounts to knowing the 
density matrix p. 

If the Q-symbol (54) determines an operator A, the values of the continuous Q- 
symbol of A at points different from those of the constellation must be functions of 
the numbers {Ai, . . . , A^J. For a coherent state |no) 7^ l^jy), not a member of the 
constellation, this dependence reads explicitly 

1 Ns 

(no|A|no) = ^^EAn(no|n.)r. (56) 

Here the P-symbol A'^""' of A is required, calculated from its Q-symbol by means of 
once the matrix 



G = G^^/= Tr ^Q^Q^>^ = \{n^\n„i)\ , (57) 

has been inverted. Furthermore, knowledge of provides immediately the dual kernel 
Q'^ via (^8]) but no explicit general expression such as ( pSf ) is known. 

It will be shown now how to directly determine the matrix elements of the dual 
kernel without using the inverse of G. The orthogonality of the kernel and its dual, Eq. 
(^^ can be written as 

1 



K = Tr 

2s + 1 



g.<5'^'J=^— Y E {m'\QM){m\Q'''\m'), (5J 

m,m'=—s 



using the completeness relation for the z eigenstates |m, n^). Introduce an (A^^ x Ng) 
matrix Q with elements Qu,mm' = ('^IQi/l^'), where the index {m,m') of the columns 
runs through A^^ values according to 

{{2s, 2s), {2s, 2s - 1), ... , (2s, 0), (2s - 1, 2s), ... , (0, 0)} . (59) 
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As is obvious from (|58|), the matrix elements of the dual kernel, Q'^^i = {m\Q'^\m') can 
be read off once the inverse of the matrix Q has been found. The expansion coefficients 
of a coherent state |n) in the z basis |m, n^) are given by 



(m,n^|n) 



1+ z 



2\s 



2s 
s — m 



1/2 



(60) 



where the complex number z is the stereographic image in the complex plane of the 
point n on the sphere. Therefore, one can write Q as a product of three matrices two 
of which are diagonal: Q = D1ND2. The diagonal matrices 



Di 
D2 



diag 
diag 



+ 

/ 2s 
I 2s — m 



1/2 



2s 
2s — m' 



1/2 



(61) 
(62) 



with p = 1, . . . , Ns, and m,m' = 0, . . . , 2s, have inverses since all diagonal entries are 
different from zero. The hard part of the inversion is due to the matrix N with elements 



N, 



2s— m / *\2s—m' 



(63) 



similar to but not identical with to the structure of a Vandermonde matrix. As discussed 
in the following chapter, particular constellations give rise to matrices N with inversion 
formulae simpler than the general one. Once N has been inverted, the matrix elements 
of the dual kernel are given by the rows of the {N s x A^g) matrix 



(64) 



For discrete Q-symbols, the kernel L in (^2]), which implements the discrete * prod- 
uct, has the form: 



-'fiuX 



Tr 



Qt^QuQx = (nA|n^)(n^|ni,)(n^|nA) 



which, by using results from can be written as 

1 



where goin^, ■ n^) 
of (|66|), one has 



42s 



1 + ■ + n,, • Ha + Ha • + in^ ■ n,, A Ux) 



2s 



\ s isA{fiuX) 



(65) 

(66) 
(67) 



1^1 + ■ n,^)/2, and, defining g{^vX) as the term in round brackets 



A(/iz/A) 



In 



' g(jwX)\ 



(68) 



Therefore, the phase A has a geometrical interpretation f^: it is the surface of the 
geodesic triangle given by the points n^, n^, nx . 
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Figure 1: Examples of nested cones, free cones, and a spiral constellation for spin 
quantum number s = 1. Each set of nine points defines an allowed constellation. 

4 Constellations 

In this section examples of specific constellations are presented for which it is possible to 
prove at least that the Gram matrix has a determinant different from zero. Furthermore, 
in some cases relatively simple expressions for the dual kernel or, equivalently, for the 
inverse of the Gram matrix G are obtained. The kernel is supposed throughout to consist 
of Ns projection operators Q,y on coherent states as given in (|53D . In other words, the 
focus is on discrete Q-symbols and the P-symbols related to them. Note that, once a 
constellation has been shown to give rise to a basis in As, the inversion of its Gram matrix 
is always possible but lengthy (already for a spin 1/2): express the matrix elements of 
in terms of the co-factors of G. Four different types of constellations will be discussed 
involving randomly chosen points, points on nested cones, on free cones, and on spirals. 



Random constellations 

As shown in |16], almost any distribution of A^^ points on the sphere gives rise to an 
allowed constellation. A random selection of directions leads with probability one to an 
invertible Gram matrix. This result shows that in an infinitesimal neighborhood of any 
forbidden constellation one can find an allowed one. 

Nested cones 

Historically, this family of constellations provided the first example of allowed constella- 
tions for both integer and half- integer spins |T^. For an integer value of s, e.g., consider 



(2s + 1) cones about one axis in space, e^, say, all with different opening angles. Dis- 
tribute (2s + 1) directions over each of these nested cones in such a way that the ensemble 
of directions on each cone is invariant under a rotation about e^ by an angle 27r/ (2s -|- 1). 
For specific opening angles of the cones, the inversion of the matrix N in (^) reduces 
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after a Fourier transformation to the inversion of (2s + 1) Vandermonde matrices of size 

(2s + 1) X (2s +1). For a half-intcgcr spin the same construction is possible except that 
the directions on different cones must also lie on different meridians. There is, in fact, 
a slight generalization of this result: the same calculation with (2s + 1) arbitrary dif- 
ferent opening angles leads to (2s + 1) generalized Vandermonde matrices with nonzero 
determinant defining thus also a allowed constellations. 

Constellations on nested cones are useful also for numerical calculations because they 
allow one to distribute Ns points in a homogeneous fashion on the surface of the sphere. 
If two points of a constellation approach each other, the determinant of the matrix G 
typically becomes very large, with a disastrous effect on numerical precision. 

Free cones 

Here is another family of constellation involving (2s + 1) cones with directions located 
on them. However, now the cones may be oriented arbitrarily (no nesting), and the 
number of directions may vary from cone to cone. For example, the number of points 

on a cone can be chosen to equal the multiplicities of the spherical harmonics Yi^ with 
/ = 2s. It is claimed that allowed constellations can be identified by taking into account 
the following properties (tested numerically for values up to s = 6): 

1. The determinant of G is zero if there are more than (4s + 1) directions on a single 



2. If there are (4s + 1) points on one cone, then another cone will contain at most 
(4s — 1) points, allowing for no more than (4s — 3) directions on the third cone, 
etc. 

3. It is necessary to have directions located on at least (2s + 1) different cones. 

For a spin 1/2 these properties will be shown to hold in the next section. The first of 
these observations can be proved for arbitrary spin s by using a particular decomposition 
of the matrix G, 



exploiting the fact that a positive definite matrix can always be written as the 'square' 
of its 'root.' A lengthy calculation involving properties of rotation matrices, Legendre 
polynomials and spherical harmonics leads to a factorization, g = dy, the first matrix 
being diagonal and having (2s + 1) different entries. 



cone. 



(69) 




Z = 0, ...,2s. 



(70) 
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each value occurring (2/ + 1) times. The second matrix has columns given by the A^^ 
lowest spherical harmonics evaluated at one of the Ng points of the constellation, 



/ yoo(ni) Foo(n2) 
ri_i(nO ri_i(n2) 



•• >00(nArJ \ 



(71) 



\Y2s2sM ^2s2s(n2) ... >^2s2s(nArJ/ 



Consequently, the Gram matrix G is invertible if and only if det y 7^ 0. The matrix ([7l|) 
can accomodate at most (4s + 1) directions on one cone, corresponding to one value of 
{} with respect to some fixed axis. The subsequent multiplicities (4s — 1), (4s — 3), . . ., 
are due to applying the same argument to the remaining subspaces with dimensions 
2(/-l) + l,2(/-2) + l,... 

In physical terms, determinant of ( |7TD is easily interpreted as a Slater determinant of 
a quantum system: it equals the (totally anti-symmetric) ground-state wave-function of 
Ng non-interacting fermions restricted to move on a sphere. The node lines of this wave 
function correspond to forbidden constellations in which the corresponding operator 
kernel is degenerate, i.e., does not give rise to a basis in As- 



Spirals 

A particularly convenient constellation is defined in the following way: let the Ng direc- 
tions be defined by Ng complex numbers points constructed out of a single point zq 
(neither of modulus one nor purely real), 

Zu = z^-\ i/ = l,...,iV,. (72) 

The points are thus located on a spiral in the complex plane. 

The matrix N defined in (|63D then turns into an (A^^ x A^^) Vandermonde matrix, 
that is, 

V.^ = xr\ iy,fi = l,...,Ng. (73) 
Its inverse is known explicitly (|74D , given, for example, in [|18|, with elements 

V;,^ = „ , 5^.-.({xA}f=i - X,) , (74) 

where SNs-ui{xx}x=i — x^) is the symmetrical function constructed out of the N — v 
numbers x^, with v ^ \. One has, for example, 5*2 (xi, X2, Xs) = X\X2 + X\X-i -|- ^2X3. 
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5 Discrete Moyal representation for a spin 1/2 



In this section the discrete Moyal representation will be worked out in detail for a spin 
with quantum number s = 1/2, allowing for explicit results throughout. For clarity, it 
is assumed from the outset that the kernel consists of four projection operators 

= |n^)(n^| , u = l,...,Ns. (75) 

It is easy to generalize the results derived below to the case of four linear combinations 
of I ± ni,)(±nj,| compatible with Eq. (^2|). 

Let us start with the determination of the dual kernel which can be found by the 
intermediate step of inverting the (4 x 4) Gram matrix with elements 

G = |K|n.)P = l(l + n^-n,) . (76) 

This matrix is easily factorized: G = g^g/2, where 



/ 


1 


1 


1 


1 


\ 




nix 


n2x 


n3x 


n4x 






niy 


n2y 


nsy 


riAy 




V 




n2z 




riAz 


) 



(77) 



The absolute value of the determinant of g is proportional to the volume of the tetra- 
hedron defined by the four points rii, on the surface of the sphere implying | det G| = 
ISVtetra- Siucc a 'flat' tetrahedron has no volume, the entire set of forbidden constella- 
tions has a simple geometric description: 

det G = <^==^ the four points n^, are located on a circle on S^. (78) 

Consequently, allowed constellations are characterized by three vectors on a cone (any 
three points on a sphere define a circle), plus any fourth vector not on this cone. This 
agrees with the earlier statements about free-cones constellations. 

Here is a simple way to invert the matrix g and subsequently G. Consider a matrix 



/ 1 


J X 




Pz\ 

Pz 


1 


r2 

J X 


r2 

Jy 


1 


fx 


Jy 


Pz 


V 1 


ft 


Jy 


n 1 



defined in terms of four vectors = (/^, fy, f^) not required to have length one. The 
matrix elements of the of product f and g are given by 

(fg);: = 1 + P ■ . (80) 
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This is a diagonal matrix if the scalar products ■ equal —1 whenever jj, ^ u. 
Geometrically, such four vectors are constructed easily: the vector points to the 
unique intersection of the three planes tangent to the sphere at the points — n2, — na 
and — n4. Analytically, this vector reads 

, n2 A ng + ng A n4 + n4 A 

I — 7 r , (oil 

(n2 A ns) ■ n4 

and the three remaining vectors follow from cyclic permutation of the numbers 1 to 4. 
With this choice the inverse of the matrix g can be written as 

g"' = d-if, (82) 

where d is the diagonal matrix in (^): (^w = 1 + f ^ ■ n,^. Consequently, the inverse of 
the Gram matrix G for a general allowed constellation is given by 

= 2d"ifftd-S (83) 

having matrix elements 

In general, the elements Q'^ of the dual kernel will thus be linear combinations of all 
four projection operators Qy. 

It is interesting to express the kernel and its dual in terms of the Pauli matrices 
(T = [ax-, (Ty, az): 

4 = ^(l + n,-a), Q^ = '^ llllZy ^^^^ 

allowing one to show easily that they satisfy the required duality. 
For reference, we give the Q- and P-symbols of the spin operator 

1 2 
= -rij, , s*^ = , (86) 



and of the identity, 

4 



" ' 1 + f ■ ' 

and the symbols of arbitrary operators for a spin 1/2 follow from linear combinations. 



7) 
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6 Discussion 



Operator kernels have been used for a systematic study of phase-space representations of 
a quantum spin s. The kernels have been derived from appropriate Stratonovich-Weyl 
postulates taking slightly different forms for continuous and discrete representations, 
respectively. Emphasis is on the discrete Moyal formalism which allows one to describe 
hermitean operators, including density matrices, by a minimal number of probabilities 
easily measured by a Stern-Gerlach apparatus. As a useful by-product a natural and 
most economic method of state reconstruction emerges when a quantum spin is described 
in terms of discrete symbols. Further, Schrodinger's equation for a spin s turns into a 
set of coupled linear differential equations for (2s + 1)^ probabilities ||19[| . 

In addition, a new form of the kernel defining continuous Wigner functions for a 
spin has been obtained p^): it has been expressed as an ensemble of operators obtained 
from all possible rotations of one fixed operator. This is entirely analogous to an elegant 
expression of the kernel for particle- Wigner functions as the ensemble of all possible 
phase-space translations of the parity operator. Therefore, continuous phase-space rep- 
resentations for both spin and particle systems now are seen to derive from structurally 
equivalent operator kernels. 

The discrete symbolic calculus is an interesting 'hybrid' between the classical and 
quantal descriptions of a spin. On the one hand, this representation is equivalent to 
standard quantum mechanics of a spin. On the other, the independent variables carry 
phase-space coordinates as labels (^,^). However, only a finite subset of points in 
phase space (corresponding to an allowed constellation) are involved reflecting thus the 
discretization characteristic of quantum mechanics. 

The Ns projections operators associated with a constellation of points define a non- 
orthogonal basis for hermitean operators acting on the Hilbert space of the spin. Each 
projection is a positive operator, and, altogether, they give rise to a resolution of unity. 
One might suspect that they define a positive operator-valued measure |^ or POVM, 
for short. However, this is not the case since the closure relation does not involve 
just the bare projections but they are multiplied with factors some of which necessar- 
ily take negative values. Such an obstruction through 'negative probabilities' is not 
surprising; other phase-space representations are based on quantum mechanical 'quasi- 
probabilities,' known to have this property, too. 

Let us close with a synopsis of the fundamental Moyal-type representations for par- 
ticle and spin systems known so far. 
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self-dual kernel 


dual pairs 


particle 


Wigner functions 
[ unknown ] 


P-, Q-symbols 
[ unknown ] 


spin 


Stratonovich/ Varilly 
[ impossible ] 


Berezin symbols 
[ P'^-, g^.-symbols] 



The table provides both a summary and points at open questions. The individual 
entries give the names of the familiar continuous phase-space representations while the 
corresponding quantities for the discrete formalism are in square brackets. Future work 
will focus on developing a discrete Moyal-typc formalism for a quantum particle. To do 
this, one must exhibit, for example, a pair of dual kernels one of which would consist 
of a countable set of projection operators on coherent states. This set is required to be 
a basis in the linear space of (bounded?) operators on the particle Hilbert-space. It is 
not obvious in which way the associate discrete P-symbol would reflect the subtleties of 
its continuous counterpart which may be singular. Similarly, the existence of a self-dual 
discrete kernel for a quantum particle is an open question. 
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